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1. Introduction 

The search for supersymmetry [1, 2] is a central part of the physics program at the proton- 
antiproton collider Tevatron with a centre-of-mass energy of = 1.96 TeV and at the 
proton-proton collider LHC, which started operation in 2010 at x/S* = 7 TeV. In particular 
squarks and gluinos, the coloured supersymmetric particles, may be produced copiously in 
hadronic collisions. The hadroproduction of top squarks (stops) [3] is an important special 
case, since the strong Yukawa coupling between top quarks, stops and Higgs fields gives 
rise to potentially large mixing effects and mass splitting [4]. The same holds, albeit to 
a lesser extent, for bottom squarks (sbottoms). Moreover, if the scalar masses in unified 
supersymmetric theories are evolved from universal values at high scales down to low scales, 
the lighter of the stop mass eigenstates is generally driven to the lowest value in the entire 
squark mass spectrum. The search for the lightest stop therefore plays a special role in the 
quest to find signals of supersymmetry at hadron colliders. 

Searches at LEP [5, 6] and the Tevatron [7] - [11] have placed lower limits on the lighter 
stop and sbottom mass eigenstates in the range between about 70-200 GeV, depending 
on the choice of supersymmetric parameters. The LHC will extend the range of sensitivity 
into the TeV-region [12, 13]. 

In the minimal supersymmetric extension of the Standard Model (MSSM) [14, 15] with 
R-parity conservation, stops are pair-produced at hadron colliders: 

pp/pp titi +X and 1212 +X , (1.1) 
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where ti and t2 denote the hghter and heavier mass eigenstate, respectively. The hadropro- 
duction of mixed tit2 or ^2^1 final states is strongly suppressed since it can only pro- 
ceed through electroweak channels or QCD-induced loop diagrams [3, 16, 17]. Sbottom 
hadroproduction is described in a completely analogous manner, so we will focus our dis- 
cussion on stops. We will, however, comment on potential differences between stop and 
sbottom hadroproduction, and provide benchmark cross sections for sbottom production 
at the Tevatron and the LHC. 

Accurate theoretical predictions for the stop-pair cross sections are crucial to derive 
exclusion limits [7] -[11] and, in the case of discovery, can be used to determine the stop 
masses and properties (see e.g. Refs. [18, 19, 20]). The cross sections for the stop-pair 
production processes (1.1) have been calculated at next-to-leading order (NLO) in su- 
persymmctric QCD (SUSY-QCD) [3]. The SUSY-QCD corrections significantly reduce 
the rcnormalization- and factorization-scale dependence and increase the cross section 
with respect to the leading-order (LO) predictions if the renormalization and factoriza- 
tion scales are chosen close to the stop mass. Electroweak corrections have been calculated 
as well [21, 22]. Although they can be sizeable at large invariant masses and large trans- 
verse momenta, they are moderate for the inclusive stop cross section. The SUSY-QCD 
calculation of Ref. [3] has been implemented in the public computer code Prospino [23] 
and presently forms the theoretical basis for the stop mass limits obtained at the Tevatron. 

A significant part of the large NLO SUSY-QCD corrections to squark hadroproduction 
can in general be attributed to the threshold region [24, 3] where the partonic centre-of- 
mass energy is close to the kinematic threshold for producing massive particles. In this 
region the NLO corrections are dominated by contributions from soft gluon emission off 
the coloured particles in the initial and final state and by Coulomb corrections due to 
the exchange of gluons between the massive particles in the final state. The soft-gluon 
corrections can be taken into account to all orders in perturbation theory by means of 
threshold resummation. A considerable amount of work has recently been devoted to 
the calculation of threshold logarithms for total gluino and squark cross sections [25] - 
[30]. Final-state stops are excluded in these calculations and all other squark flavours, the 
so-called light-flavour squarks, are treated as being mass degenerate, neglecting possible 
mixing effects. 

In this work, we extend the previous analyses of threshold resummation for the hadro- 
production of gluinos and mass-degenerate light-flavour squarks at next-to-leading loga- 
rithmic (NLL) accuracy [25, 27, 30]. Firstly, we consider the hadroproduction of stops 
and non-mass-degenerate sbottoms. Secondly, we study the impact of NLO and NLL 

corrections on the transverse-momentum distributions. Since theoretical predictions for 
differential distributions are input to the experimental analyses, it is important to assess 
how the shape of the distributions is affected by higher-order corrections. The thresh- 
old resummation for transverse-momentum distributions has been studied extensively for 
Standard- Model processes, see e.g. Refs. [31] -[37], but not yet for SUSY processes. 

The paper is structured as follows. In section 2 we review the basic features of stop and 
sbottom hadroproduction, and we briefly discuss the application of threshold resummation 
to the transverse-momentum distribution. Section 2 also contains the specific stop-pair 
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formulae that enter the calculation of the resummed cross sections. State-of-the-art SUSY- 
QCD predictions for stop hadroproduction at the Tevatron and the LHC, including NLO 
corrections and NLL threshold resummation, are presented in section 3. We discuss the 
impact of the NLO-I-NLL corrections on total cross sections and transverse-momentum 
distributions and provide an estimate of the theoretical uncertainty due to scale variation 
and the parton distribution functions. We conclude in section 4. The dependence of the 
stop and sbottom cross sections on the choice of supersymmetric parameters can be found 
in the appendix, where also some predictions for specific benchmark scenarios are given. 

2. Stop and sbottom pair production 

Let us first review some basic features of the stop and sbottom pair-production cross 
sections. At LO the hadroproduction of stop pairs proceeds through quark-antiquark 
annihilation and gluon-gluon fusion: 

qq — )■ titi and ^2^2 , 

gg — titi and t2t2 ■ (2-1) 
The corresponding Feynman diagrams are shown in Fig. 1. In contrast to the hadropro- 




Figure 1: Leading-order Feynman diagrams for the production of stop pairs through quark- 
antiquark annihilation (first line) and gluon-gluon fusion (second line). 

duction of light-flavour squarks, no t-channel gluino-exchange graph occurs in the quark- 
antiquark channel. In such a t-channel graph the initial-state quarks should have the same 
flavour as the final-state squarks, but since top quarks are excluded as initial-state partons, 
the gluino-exchange graph is absent. 

The mass eigenstates ti and t2 are related to the weak interaction eigenstates tL and 
ifi through mixing: ii = ii cos % -|- sin % and ^2 = —tL sin % -|- iji cos %. The masses 
rrij^, rrij^ and the mixing angle Oj; are obtained from diagonalizing the stop mass matrix 
and are determined by Standard-Model and soft-supersymmetry-breaking parameters [4]. 
As mentioned in the introduction, mixed pairs tit2 or t2ti cannot be produced in lowest 
order QCD since the gtt and ggtt vertices are diagonal in the chiral and in the mass basis. 
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The description of sbottom hadroproduction pp/pp — > bibi and 62^2 is completely 
analogous to that of stop production. The only differences occur in the bb — > bb channel, 
where the initial-state bottom quarks do allow a t-channel gluino-exchange graph that 
gives rise to extra contributions. As will be discussed below, these contributions lead to a 
less-suppressed LO threshold behaviour than the s-channel gluon-exchange contributions. 
However, we shall demonstrate in the appendix that their numerical impact on the hadronic 
cross sections is negligible. Thus, for all practical purposes, the LO and higher-order cross- 
section predictions obtained for stop-pair production apply also to sbottom-pair production 
if the input parameters, i.e. masses and mixing angles, are modified accordingly. 

When decomposed into s-channel colour states, the LO partonic cross sections for the 
subprocesses (2.1) read in generic notation: 

C7^°^ -= = , (2.2) 

(0) _ a^TTiN^ - 1) 3 .2 3) 



c 



(0) "s^ 

a — 



gg^ii,! " Nc{N^ - l)s [\2^ ^ ^ \s vTT^ 
gg^tt,8 2 gg^iiA 

Sm? \ ( m? 2m^ \ / 1 — /? 



(2.4) 



( 1 8m \ ^ ( m 2m \ , 
(l2 + ^)^ + (v + ^)^°^ 



1 + 



(2.5) 



where ctg is the strong coupling, ,s the invariant partonic energy and A^^ the number of 
colours. The colour labels 1 and 8 refer to the familiar singlet and octet colour states 
in SU(3), but all analytic results in this paper are derived for a general SU(iVc)-theory 
in order to make the colour structure more transparent. Also note that we use a generic 
notation for the final-state particles and the associated kinematics: i can be either of the 
two stop mass eigenstates, with m being the corresponding mass and P = a/1 — Am^ / s 
the corresponding velocity. 

The expressions (2.4) and (2.5) for the gluon-gluon fusion cross section agree with 
the corresponding expression for light-flavour squark production, cf. Ref. [24]. However, 
the quark-antiquark annihilation contribution (2.3) is different due to the absence of the 
t-channel gluino-exchange graph. As a consequence, the LO cross section for qq — > it 
proceeds through s-channel gluon exchange only and is proportional to , as opposed to ^ 
for other squark flavours. This behaviour is the combined effect of the standard phase- 
space suppression factor /3 and an additional P-wave suppression oc near threshold: the 
pair of scalar particles needs to be produced in a P-wave state to balance the spin of the 
intermediate gluon. 

Note that the LO cross sections (2.2)-(2.5) only depend on the mass m of the produced 

stops and not on any other supersymmetric parameters. At NLO, however, the stop mixing 
angle 9i enters through corrections involving the ttg vertex and the four-squark couplings. 
As a result, already the analytical expressions for the iiii and ^2*2 NLO cross sections 
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are different. Furthermore, virtual corrections involving squark and gluino loops introduce 
a dependence of the stop-pair cross section on the masses of the squarks and the gluino. 
The dependence on all these other supersymmetric parameters turns out to be mild, as 
illustrated in the appendix. 

In the threshold region the NLO stop-pair cross section is dominated by soft gluon 
emission, which leads to corrections oc log* /3 (i = 1,2), and Coulomb corrections oc 
In the notation of Ref. [3] the threshold behaviour of the total NLO cross sections reads^: 



a 



(l,thr) 
qq—^tt 



+ 



48m2 
2Cf 



47r2 



log2(8^2) _ ^ iog(8^2) _ i^g(8^2) log(4) 
6 Vm^/ 



2Cf — Ca 



(2.6) 



and 



a 



(l,thr) ^ 
gg-^ti 



+ 



16m2 

2Ca 
47r2 



■/3 l + 47ras(/x2 



16/3 



log2(8/32) - 41og(8/32) - log(8/32) \og(^ 



- 2 47r2 



log(8;32) 



(2.7) 



Here n denotes the factorization and renormalization scales, which we keep equal in this 
analysis. The first line in both expressions contains the gluon corrections that can be 
attributed to the final-state particles, while the terms in the second line, proportional to 
either Ci;' = (A^^ — l)/(2A'^c) or = A^c correspond to initial-state radiation from a quark 
or a gluon line, respectively. The final-state radiation consists of two parts. The logarithmic 
soft-emission terms are proportional to the colour charge of the final state [38, 25, 28, 30] 
and are therefore absent for the singlet colour state and proportional to Ca for the octet 
colour state. As a result, in the expression for the gluon-gluon fusion channel these terms 
are multiplied by (A"^ — 4)/{N^ — 2), which is the ratio of the colour-octet and the total 
cross section at threshold (cf. [38]). The Coulomb corrections receive contributions with 
colour factors Cp and the total colour charge of the final state. Note, finally, that the 
threshold behaviour of the cross section is determined by QCD dynamics and does not 
involve any supersymmetric parameters other than the mass of the produced stops. 

2.1 Threshold resummation for the total inclusive cross section 

In this paper we shall improve the NLO prediction for the hadroproduction of stop and 
sbottom pairs [3] by resumming the leading and next-to-leading threshold logarithms of 
the form 

a^log"/?^ , m<2n (2.8) 



^Ref. [3] contains a misprint in the log(8/3^)-coefRcient for the gg-channel. In Eq. (9) of [3] the coefHcient 
-107/(367r^) should be replaced by -155/(367r^). 
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to all orders, and to next-to-lcading logarithmic (NLL) accuracy. The resummation forma- 
lism developed for heavy-quark production [39, 38, 40] can be applied directly to squark- 
antisquark production as the colour structure is identical for both cases. The resummation 
is performed in the space of Mellin moments (N-space) following the procedure outlined in 
Refs. [25, 27, 30], where the hadroproduction of light-flavour squarks and gluinos at NLL 
accuracy has been considered. 

From the viewpoint of threshold resummation there is a notable difference between 
the case of stops and light-flavour squarks, arising from the extra suppression of the 
quark-antiquark annihilation cross section cr^^^ii near threshold. In AT-space this effec- 
tively produces an extra factor \/N compared to the light-flavour squark case, resulting in 
leading contributions oc log* (iV)/iV instead of log' (AT). For the seemingly analogous case 
of threshold resummation for the deep-inelastic structure function Fl, which has a leading 
behaviour of the type log*(A^)/A^ as well, differences from the expected NLL resummation 
structure were revealed in Refs. [41] and [42]. However, the NLL resummation procedure 
developed for light-flavour squarks applies to stop production as well. In the case of Fl the 
extra factor 1/iV with respect to the N dependence of the structure function F2 arises due 
to a special initial-state jet function [43, 44] associated with the longitudinal projection. 
In contrast, in our case the log*(A^)/A^ dependence arises from projecting onto the P-wave 
final state, which docs not depend on the initial state jet function. Moreover, at ©(cts) lead- 
ing and subleading \og{N)/N corrections can be computed from the one-loop calculations, 
and they do exhibit the pattern expected from NLL threshold resummation. Finally, in 
view of the different threshold behaviour of the qq and gg channels one might worry about 
the possibility that these channels mix in the NLL threshold resummation as a result of 
soft-quark emissions, which is not the case for top-quark and light-flavour squark produc- 
tion. However, we have checked that to NLL accuracy no such mixing occurs at 0{I3^). 
Based on these observations we are confident that the expressions of Refs. [25, 27, 30] can 
be applied to inclusive stop-pair production as well. 

The new elements that enter the NLL resummed cross section for stop pair production 
are the LO partonic cross sections in iV-space (indicated by a tilde), decomposed into 
s-channel singlet and octet colour states. They are given by: 

-= ( AT) = , (2.9) 

qq-^tt,l^ ' ' ^ ' 

a^°> - (N) = —5 - '-^'^ ^ ('2 10) 

"qq-^ilsS'^' 64m2 r(Ar + 7/2) ' ^ ' 

^(0) = 1 NiN + 3) + 4 TjN + l) 

°gg-^tt,i^ > 16m2 7V^(iV2-l) (Ar + 2)(A^ + 3) r(Ar + 5/2) ' ^' ' 

"gg^ttfi^ > 2 64m2 A^2_i (7v + 2)(A^ + 3) r(iV + 7/2) ^ ' 

The results for the gluon-gluon channels (2.11) and (2.12) agree with those presented in 
[27] for light-flavour squarks and = 3, and are reproduced here for completeness for 
general SU(iVc). Note that the results in [27] include an additional factor 2nf from the 
summation over the L and R squarks of nj light flavours. 



-6- 



2.2 Threshold resummation for the transverse-momentum distribution 

Similar to the inclusive cross section, soft-gluon corrections can also dominate the transverse- 
momentum distribution of the stops if the production takes place sufficiently close to the 
edge of phase space. We now briefly review the general construction of threshold-resummed 
transverse-momentum distributions (cf. [31]- [37]) and its application to stop-pair produc- 
tion. 

We start with the general framework applicable to the hadroproduction of a pair 
of massive coloured particles. Assuming that one of the final-state particles carries a 
transverse momentum pT, the minimal energy necessary to produce the system is 2mT, 
where the transverse mass rriT is defined by 



rriT 



\lm?+pl . (2.13) 



For the py-distribution the dominant contributions originating from soft gluon emission 
have again the structure (2.8), with the variable ^ replaced by (cf. [45]) 

I5t = ^Jl-iml/s . (2.14) 

Like in the case of the total inclusive cross section, the resummation of the logarith- 
mic threshold corrections to pf-distributions takes place in the space of Mellin moments. 
However, in this case the Mellin transform is taken with respect to the scaling variable 
x"^ = Am^/s, i.e. 



da. 



dpT 



{N,pT,n') = rd4(4)^-^^^(4,PT,M') (2.15) 

Jo (f'PT 



for a generic partonic subprocess ij — > kl. 

In A^-moment space the resummcd partonic pr-distribution for the hadroproduction 
of two massive coloured particles is given by 

-^i^^PT^^') = E -^^{N,PT,^')Qj^klAN,PT,^') (2.16) 

X MN+l, m') A,(A^+1, Q^ //') A^_^,i j{N + 1, Q^, m") , 

with Q = 2mT representing the resummation scale and the index / indicating all possible 
colour states of the hard scattering. To NLL accuracy, the coefficients Cij^kij , Pt , 1^"^) 
are equal to 1 for all channels and colour structures. The functions Aj and Aj sum the 
effects of the (soft-)collinear radiation from the incoming partons. They are process- and 
colour-independent, and are therefore a universal ingredient in all threshold-resummed 
cross sections. Explicit expressions for the Aj can for instance be found in Ref. [27]. The 
function A^^^^^ j accounts for soft wide-angle emission and depends on the production 
process and the colour channel. At NLL accuracy it is given by 

AgL,.(iV, Q^ ^^') = exp [ f''^'' ^ ^ A,^.,. ] ■ (2.17) 
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The values of the coefficients Dij^^ij follow from the threshold limit of the one-loop 
soft anomalous-dimension matrix in the way described e.g. in Ref. [30]. If an s-channel 
colour basis is used, the soft anomalous-dimension matrix for the p-r-distribution becomes 
diagonal in colour space in the threshold limit ^ 0, leading to the colour-diagonal form 
of the resummation formula given in Eq. (2.16). This is similar to threshold resummation 
for the total cross section, where the soft anomalous-dimension matrix becomes diagonal 
in colour space in the corresponding threshold limit /3 — t- 0. Since the threshold limit is 
defined differently for the total inclusive cross section and the pT-distribution, the values 
of the D-coefficients are different as well; in particular, the Z)-coefficients now carry pr- 
dependence. 

For the stop-pair production processes we have 

D,,/,,-.m = -^C^{^eL^^ + l) , (2.18) 



r 2 

^,,M-.«,8 = - (ReL^, + 1)+Ca log(^) + ReL^, 



(2.19) 



with 



and 



n.L. = l±Mlo,C^) (2.20) 



2/3r v1 + /3t 



Pt = lim Pt = Jl-rri^/ml . (2.21) 

These expressions agree with those obtained for heavy-quark production in the context of 
joint resummation [37]. 

In addition to the soft radiative factor, the other new elements which have to be calcu- 
lated in order to obtain resummed predictions from Eq. (2.16) are the colour-decomposed 
LO pr-distributions in AT-moment space. They are obtained from 



(2.22) 
(2.23) 
(2.24) 



qq—>tt,l 

dpT 


= , 






da^'^ - 

qq-^tt,S 


2 


2(iV2 - 1) p3 




dpT 


S3 






da^'^ 


2 


2 PT 


4 1 4 

m + pj. 


dpT 




N,{N^ - 1) Pt 


4 ' 


da^'^ -= 
gg-^tt,8 


Ar2 - 


- 4 gg-^tt,l 1 


OgTr Nc 


dpT 


2 


dpT 


S2 iV2- 



4 I 4 

■PtPt 2.25 



Taking a Mellin transform of (2.22)-(2.25) with respect to the threshold variable one 
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finds 



da 



da 



(0) 

qq—^tt,l 

dpT 

(0) 



{N) = 0, 



dpi 



{N) = 



a 



f 7r3/2 N^-l TjN + 3) ^ 
r{N + 7/2) m| ' 



32 



T{N + 2) prim^ + p^) 



da 



dpi 

(0) 

dpT 



(N) 



8 Nc{N^ -1) r{N + 5/2) 

t\t2 a da^^^ -.= ^2„3/2 



r(iV + 2) pT(m^+pf) 



dpT 



32 Ar2 _ 1 T{N + 7/2) 



(2.26) 
(2.27) 
(2.28) 

(2.29) 



Having calculated the resummed partonic expression in A^-space, Eq. (2.16), the re- 
summed hadronic pr-distribution is obtained by the inverse Mellin transform 



da 



(res) 



dpi 



{xt,pt,h^) = 



^,3=1,(1,9 



dN 



(4)"%/^l(^^+l,M')/,7/..(^^+l,M') 



,~(rcs) 



(2.30) 



where x'^ = Am^/ S is the hadronic scaling variable. The functions fijh^ and fj/h2 ^'^^ 
Mellin moments of the parton distribution functions for the initial-state hadrons h\ and 
h2- In order to retain the information contained in the NLO pr-distributions [3], the NLO 
and NLL results are combined through a matching procedure that avoids double counting 
of the logarithmic terms in the following way: 

(NLO+NLL) (NLO) 
hih^-^kl 12 2\ "-"hihi-^kr 12 2\ 



dal 



dpi 



dpi 



'-,J=<1,<1,9 



^^(res) 



dpT 



{N,pT,fi^) - 



da, 



(res) 

ij—^kl 



{N,pT,l?) 



(NLO) 



(2.31) 



where (da^^^^) /dpr) | (NLO) represents the perturbative expansion of the NLL pT-distribution 
(2.16) truncated at the order of associated with NLO. To evaluate the inverse Mellin 
transform in Eqs. (2.30) and (2.31) we adopt the "minimal prescription" of Ref. [46] for 
the integration contour CT. 



3. Numerical results 

We present numerical results for the NLL-resummcd cross sections and transverse-momen- 
tum distributions matched with the complete NLO results for stop-pair production at the 
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Tevatron (-v/S* = 1.96 TeV) and the LHC {y/S = 7 and 14 TeV). The matching procedures 
for the total cross section and for the transverse-momentum distribution are described in 
Ref. [30] and in Eq. (2.31), respectively. We use the notation NLO+NLL for matched 
quantities in the following. The NLO corrections are calculated using Prospino [23], based 
on the calculation presented in Ref. [3] . The QCD coupling ag and the parton distribution 
functions (pdfs) at NLO are defined in the MS scheme with five active flavours. The mass 
of the stop is renormahzed in the on-shcll scheme and the SUSY particles are decoupled 
from the running of ctg and the pdfs. Since mixing enters explicitly only through higher- 
order diagrams, the angle 0^ need not be renormahzed and one can use the lowest-order 
expression derived from the stop mass matrix. 

As our default, hadronic total cross sections and transverse-momentum distributions 
are obtained with the 2008 NLO MSTW pdfs [47] and the corresponding a^iMz) = 0.120. 
The NLL corrections are convoluted with pdfs in Mellin space, derived with the program 
PEGASUS [49] based on the MSTW parametrization at the initial factorization scale. For 
the total hadronic cross sections, we have also used the method of Ref. [50] to evaluate 
the NLL cross section with standard parametrizations of pdfs in x-space. We find however 
a much better numerical stability when the total NLL cross sections are evaluated with 
Mellin-space pdfs. 

Beyond LO the cross section does not only depend on the stop mass, but also on the 
gluino mass rrig, the average mass of the first and second generation squarks niq and the 
mixing angle 9^. For this reason we have adopted the SPSla' benchmark scenario [51] for 
our numerical analysis. Taking a top-quark mass of mt = 172.5 GeV [52] and as{Mz) = 
0.120 the SPSla' scenario corresponds to m~g = 610 GeV, rrig = 560 GeV, sin(26'j) = 0.932 
and stop masses of = 367 GeV and m^^ = 590 GeV [53]. However, in order to focus 
on the mass dependence of the cross section and the NLO+NLL corrections, we vary the 
mass of the stop while keeping the other SUSY parameters fixed. As shown in Ref. [3] 
and discussed in more detail in the appendix, the dependence of the cross section on the 
additional SUSY parameters is small, justifying this procedure. Note that the numerical 
results presented for stop production also apply to sbottom production when the same 
input parameters are adopted. In the appendix we show that the impact of bottom-quark 
induced contributions to sbottom hadroproduction is negligible and present benchmark 
predictions for the sbottom cross section. 

3.1 Results for the total cross section 

Let us first discuss the scale dependence of the SUSY-QCD total cross-section prediction. 
Fig. 2 shows the scale dependence in LO, NLO and NLO+NLL for titi production, using 
rrif^ = 200 GeV and 500 GeV at the Tevatron and LHC, respectively. Here and in the 
following, we present results for the LHC operating at 7 TeV and at 14 TeV centre-of-mass 
energy. Note that the LO predictions are obtained with LO pdfs and the corresponding LO 
values for as [47]. The renormalization and factorization scales arc identified and varied 
around the central scale fio = m^^ from fi = fio/W up to fi = bfiQ. We observe the usual 
strong reduction of the scale dependence when going from LO to NLO. A further significant 
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improvement is obtained when the resummation of threshold logarithms is included, in 
particular for stop production at the Tevatron and at the LHC running at 7 TeV. 

Near the central scale fi = rrif^ the cross section is enhanced by the SUSY-QCD 
corrections at NLO and NLO+NLL. The size of the i^-factors K^lo = cnlo/clo and 
-^NLL = ctnlo+nll/cnlo Strongly depends on the stop mass and the collider, as is shown 
in Fig. 3. At the Tevatron, where the cross section is dominated by gg-annihilation for 
large stop masses, the NLO X-factor is moderate and ranges from roughly 1.2 to 1.03 
for stop masses in the range between 100 and 300 GeV. A further enhancement by up to 
7% is found for large stop masses when the NLL resummation is included. At the LHC, 
the gg initial state is dominant and the QCD corrections are in general larger. For iiii 
and ^2^2 production at the LHC we consider the mass ranges 100 GcV < rn^^ < 550 GeV 
(lower horizontal axis) and 550 GcV < m^^ < 1 TeV (upper horizontal axis). At 7 TeV we 
find NLO corrections ranging from about 40% at the lower end of the stop mass range to 
about 20% for stop masses near 1 TeV. The NLL resummation leads to a further increase 
of the cross-section prediction of approximately 10% for stop masses in the TeV-range. At 
14 TeV centre-of-mass energy the NLO corrections to stop production are significant and 
increase the LO cross section by around 35% for moderate ti masses and by up to 40% 
for t2 with nip^ ~ 600 GeV, while the impact of the NLL resummation is modest with 
at most 5% further increase for stop masses in the TeV-rangc. The singularities at the 
stop-decay threshold rrif = nit + iT^g = 782.5 GeV originate from the stop wave- function 
renormalization. They are an unphysical artefact of the on-shell approach of Ref. [3] and 
could be removed by taking into account the finite widths of the unstable stops. Note that 
the NLO cross sections for titi and ^2^2 production are not identical, even if the masses 
are taken equal. The reason for this is that the stop mixing angle contributes in different 
ways to both reactions at NLO, as discussed in section 2. Furthermore, while we vary 
the mass of the stop particle that appears in the final state, the mass of the other stop, 
which enters the loop corrections, is set to its SPSla' value and thus differs for ii and t2- 
However, numerically the difference between the two NLO total cross sections is moderate. 
The NLL resummation does not involve any SUSY parameters apart from the stop mass 
itself and thus affects the titi and t2t2 cross sections in the same way. The NLL ii'-factors 
have a tiny SUSY-parameter dependence, which enters through the ratio cnlo+nll/cnlo- 

Predictions for the LO, NLO, and NLO-I-NLL total cross sections are shown in Fig. 4 
and Tables 1-3 for iiii production at the Tevatron and titi 7*2*2 production at the LHC 
with 7 TeV and 14 TeV centre-of-mass energy. In fact, the cross sections for iiii and ^2*2 
production at equal masses are indistinguishable on the scale of Fig. 4. We thus refrain 
from showing additional plots for *2*2 production. The results shown in Fig. 4 represent the 
state-of-the-art SUSY-QCD predictions at NLO+NLL accuracy. The error bands include 
the NLO-I-NLL scale variation in the range < /x < 2m£^ as well as the NLO pdf 

uncertainty, added in quadrature. The pdf uncertainty is obtained with the help of the 
90% C.L. MSTW error pdfs [48]. More detailed information is available in Tables 1-3. 

In Table 1 we present results for titi production at the Tevatron. As discussed before, 
we observe an increase of the cross-section prediction near the central scale when going from 
LO to NLO and a further enhancement when NLL threshold resummation is included. The 
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scale dependence in the range < n < 2m^^ is reduced from about ±50% at LO to 

about ±10% at NLO±NLL. The estimated pdf uncertainty is approximately 5%. We also 
present cross-section predictions obtained with the CTEQ6 pdf set [54] and an estimate 
of the corresponding pdf error. The difference between the MSTW and CTEQ results is 
particularly pronounced at large stop masses, m^^ 300 GeV, where the cross sections 
obtained with CTEQ pdfs are about 7% larger than the ones obtained with MSTW pdfs. 
We observe that the CTEQ pdf error estimate of about 10% is roughly twice as large as 
that of MSTW. 

Results for iiii and t2i2 production at the LHC with 7 TeV are collected in Table 2. 
Here the impact of the SUSY-QCD corrections at NLO is large, with the NLL resummation 
adding a further enhancement of up to 10%. The scale uncertainty of the NLO+NLL 
prediction is reduced to a level of about 10%. Unfortunately the pdf error is sizeable, 
in particular at stop masses in the TeV-region, where we find a pdf error of about 20% 
from MSTW. Also the difference between MSTW and CTEQ is significant for large stop 
masses, with a 25% increase in the prediction for stop masses near 1 TeV when going from 
MSTW to CTEQ pdfs. As before, we find a pdf error from the CTEQ analysis that is 
about twice as large as that of MSTW and reaches 45% for ~ 1 TeV. Of course, the 
large pdf uncertainty is not a specific feature of stop production. It rather generically 
affects predictions for TeV-scalc particle production at the LHC with 7 TeV ccntrc-of-mass 
energy, since these predictions are particularly sensitive to the gluon pdf at large x (see 
e.g. Ref. [55]). The conclusion therefore is that more accurate pdf determinations are 
needed in order to allow for a precise prediction of heavy-particle production during the 
initial phase of the LHC at 7 TeV centre-of-mass energy. 

Going from 7 TeV to 14 TeV at the LHC, we observe in Table 3 a significant increase 
in the predicted cross section of about a factor of 4 for stop masses around 100 GeV and 
up to a factor of about 60 for masses in the TeV-region and MSTW pdfs. Just like at 
7 TeV, the scale uncertainty of the NLO+NLL prediction is down to a level of about 10%. 
The pdf uncertainty is more moderate than at 7 TeV, ranging from 3% at small masses 
to about 10% at large stop masses for MSTW, and correspondingly from 3% to 20% for 
CTEQ. 

3.2 Results for the transverse-momentum distribution 

Let us now turn to the discussion of the transverse-momentum distributions. Figure 5 
shows a comparison between the LO, NLO and NLO+NLL distributions normalized to 
unity. We use normalized distributions in order to be able to directly read off the NLO- 
and NLL-induced changes in the shape of the distribution. As for the previous results, we 
have used the stop mass as the central scale, = m, in Figure 5. This is a possible choice as 
we do not consider regions where pr ^ m and where a p^-dependent scale would have been 
mandatory. As already observed in Ref. [56], the transverse momentum carried away by 
hard gluon radiation in higher orders softens the NLO transverse-momentum distribution 
with respect to the LO distribution. This effect is particularly visible at the Tevatron and 
at the LHC with 7 TeV centre-of-mass energy. The NLL soft-gluon resummation, on the 
other hand, does not affect the shape of the distribution significantly. To elucidate the 



- 12 - 



impact of the higher-order corrections more clearly, we display the transverse-momentum 
dependence of the NLO and NLL X-factors in Figure 6, this time using the transverse 
mass rriT = \J fri^ + Pj^ as the scale. The significant softening of the transverse-momentum 
distribution at NLO at the Tevatron and the 7 TeV LHC is reflected in the variation of the 
if-f actor, with -fCNLO dropping from roughly 1.8 at small pT to a value near one at pT ~ 2m. 
In comparison, the impact of the NLL resummation is small. A similar behaviour, albeit 
less pronounced, is observed at the LHC with 14 TcV ccntrc-of-mass energy. It would 
be interesting to see if using NLO+NLL transverse-momentum distributions would affect 
the experimental analyses, which so far have been based on LO Monte Carlo predictions. 
In this context we recall that the shape of the stop rapidity distribution is not changed 
significantly by higher-order corrections, see Ref. [56]. 

4. Conclusions 

In this paper we have performed the NLL threshold resummation for stop and sbot- 
tom hadroproduction, considering both the inclusive cross sections and the transverse- 
momentum distributions. As the lighter stop and sbottom mass eigenstates are generally 
predicted to be the lightest strongly interacting SUSY particles, the search for these parti- 
cles plays a special role in the quest to find signals of supersymmetry at hadron colliders. 

Results have been given for the Tevatron and for the LHC running at both 7 TeV and 
14 TeV centre-of-mass energy. Compared to the NLO predictions for the total cross section, 
the NLL corrections lead to a significant reduction of the scale dependence and increase 
the cross section by up to 10% for masses in the TeV range if the renormalization and 
factorization scales arc chosen close to the mass of the final-state particles. We have also 
studied the SUSY parameter dependence of the stop and sbottom cross sections and find 
small variations of at most 2%. The size of bottom-induced contributions to sbottom pair 
production is negligible numerically so that predictions obtained for stop-pair production 
also apply to sbottom-pair production when the same input parameters are adopted. 

Since px cuts are used extensively in experimental analyses, which at present are based 
on LO Monte Carlo simulations, it is important to investigate how the NLO-I-NLL matched 
corrections affect the transverse-momentum distributions. We find that the NLO+NLL 
corrections can change the shape of the pr distribution considerably and thus generally 
cannot be taken into account by using a simple K-factor. 

The NLO-I-NLL matched cross sections and pT distributions presented in this paper 
constitute the state-of-the-art QCD predictions for stop and sbottom production and can 
be employed to improve current and future searches at the Tevatron and LHC. 



-13- 



Acknowledgment s 



We would like to thank Sven Moch, James Stirling, Robert Thorne and Andreas Vogt for 
valuable discussions and the MSTW collaboration for providing us with parton distribution 
functions with improved accuracy in the evolution at large x. This work has been supported 
in part by the Helmholtz Alliance "Physics at the Terascalc" , the DFG GraduicrtenkoUeg 
"Elementary Particle Physics at the TeV Scale" , the Foundation for Fundamental Research 
of Matter (FOM), the National Organization for Scientific Research (NWO), the DFG 
SFB/TR9 "Computational Particle Physics", and the European Community's Marie-Curie 
Research Training Network under contract MRTN-CT-2006-035505 "Tools and Precision 
Calculations for Physics Discoveries at Colliders". SB would like to thank the Nikhef 
Theory Group, and MK would like to thank the Institute for High Energy Physics (IFAE) 
at the Universitat Autonoma de Barcelona and the CERN TH unit for their hospitality. 
AK acknowledges the hospitality of the Institute of Theoretical Physics at the University 
of Warsaw. 



-14- 



0.60 



0.50 



0.40 



0.30 



0.20 



0.10 



0.00 



0.1 



a i^pp ^ hti + X j [pb] 

v/S = 1.96TeV 

Ho = mt\ = 200 GeV 




NLO + NLL 
NLO 
LO 



0.2 



0..5 1 



0.100 
0.090 
0.080 
0.070 
0.060 
0.050 
0.040 
0.030 
0.020 
0.010 



0.000 



0.1 



VS = 7TeV 

Ho = vif^ = 500 GeV 



NLO + NLL 

NLO 

LO 



0.2 



0.5 1 



1.20 
1.10 
1.00 
0.90 
0.80 
0.70 
0.60 
0.50 
0.40 
0.30 
0.20 
0.10 



0.1 



f7(^pp^flfl+Xj[pb] 

x/S = 14TeV 

Ho = rrij- = 500 GeV 



NLO- 

NLO 

LO 



-NLL 



0.2 



0.5 1 



Figure 2: Scale dependence of the LO, NLO and NLO+NLL cross section for stop-antistop pro- 
duction at the Tevatron and the LHC. 
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Figure 3: The NLO and NLL if-factors for stop-antistop production at tire Tevatron and tire LHC 
as a function of tire stop mass. The scale has been set to the stop mass, i.e. fi — . 
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Figure 4: Total NLO+NLL stop-pair cross section at tire Tevatron and the LHC as a function 
of tire stop mass. The error band corresponds to the scale and pdf uncertainty of the prediction, 
added in quadrature. 
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Table 1: The LO, NLO and NLO+NLL cross sections for stop-antistop production at the Tevatron 

(\/5=1.96 TeV), including errors due to scale variation (Acr^) in the range mi 12 < fi < 2m, j: . 
Results are shown for two pdf parametrizations (MSTW08 and CTEQ6) with the corresponding 
90% C.L. pdf error estimates. 
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Table 2: The LO, NLO and NLO+NLL cross sections for stop-antistop production at the LHC 

(\/S'=7 TcV), including errors due to scale variation (Acr^) in the range 771^/2 < /i < 2toj. Results 
are shown for two pdf parametrizations (MSTW08 and CTEQ6) with the corresponding 90% C.L. 
pdf error estimates. 
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Table 3: The LO, NLO and NLO+NLL cross sections for stop-antistop production at the LHC 

(\/S'=14 TcV), including errors due to scale variation (Acp) in the range TOj/2 < p < 2mf. Results 
are shown for two pdf parametrizations (MSTW08 and CTEQ6) with the corresponding 90% C.L. 
pdf error estimates. 
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Figure 5: Normalized LO, NLO and NLO+NLL transverse-momentum distributions for stop- 
antistop production at the Tevatron and the LHC for jj, = mi^ . 
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Figure 6: Transverse-momentum dependence of the NLO and NLL ii'-factors for stop-antistop 
production at ttie Tevatron and the LHC for /x = nirp f^ . 
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A. SUSY parameter dependence of stop and sbottom cross sections 

In this appendix we shall investigate the dependence of the NLO+NLL stop and sbottom 
cross-section predictions on the supersymmetric model parameters that enter beyond LO, 
i.e. the mixing angle and the masses of the light-flavour squarks and the gluino. In the 
case of sbottom production we shall in addition quantify the impact of the bottom-quark- 
induced reaction channels bb — )• b^bk and bb/bb — )• b^bk/b^bk, involving contributions from 
t-channel gluino exchange. It will be demonstrated that the contributions of these partonic 
reaction channels are strongly suppressed numerically. Thus, for all practical purposes, 
cross-section predictions obtained for stop-pair production also apply to sbottom-pair pro- 
duction when the same input parameters are adopted. 

As in the main body of the paper we choose the SPSla' benchmark scenario [51] as our 
default. The SPSla' masses and mixings relevant for stop and sbottom hadroproduction are 
collected in Table 4. Note that the stop and sbottom masses predicted by the SPSla' and 



SPSla^ 





367 GeV 




509 GeV 


rriq 


560 GeV 




590 GeV 




549 GeV 


TJlg 


610 GeV 


sin2% 


0.932 


sin 26*^ 


0.652 







Table 4: Masses and mixings for the SPSla' benchmark scenario [51] obtained using SPheno [53] 
with the Standard Model input parameters mt = 172.5 GeV and as{Mz) = 0.120. 

other commonly used benchmark scenarios are beyond the reach of the Tevatron searches, 
as the corresponding production cross sections are too small. The SPSla' NLO-t-NLL cross 
sections for stop and sbottom production at the LHC are collected in Table 5. 
Prom the cross-section predictions one can 
conclude that only the lighter of the SPSla' 
stop mass eigenstates might be detected 
during the initial phase of LHC data taking 
at 7 TeV with 1 fb~^ of integrated lumino- 
sity. Also a dedicated search for sbottom 
production in SPSla'-like scenarios will on- 
ly be possible with higher LHC energies. 

Therefore, to address the SUSY-para- 
meter dependence and to study the impact 
of bottom-quark-induced sbottom-pair pro- 
duction, we consider two different scena- 
rios that are within the reach of the Teva- 
tron and the early LHC phase. As we did 
in the main body of the paper, we use stop and sbottom masses of 100 (200) GeV at the 
Tevatron and 100 (400) GeV at the LHC, respectively, and present results for various choices 
of the mixing angle and the light-flavour squark and gluino masses, see Table 6. Note that 
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Table 5: NLO-hNLL SUSY-QCD cross sections 
for stop and sbottom pair production at the LHC 
for the SPSla' benchmark scenario. The MSTW 
pdfs have been adopted and the scale has been set 
to the mass of the particles produced. 



-23- 



the NLL resummation only depends on the final-state particle mass. The dependence on 
the other SUSY parameters enters exclusively through the NLO virtual corrections. The 
numbers hsted in Table 6 reveal that the dependence of the cross section on the mixing 
angle, the gluino mass, and the light-flavour squark masses is small indeed, with variations 
of at most 2% both at the Tevatron and at the LHC. 
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15.9 


0.255 


425 


0.222 


1.78 X 10^ 


2.39 





15.9 


0.254 


423 


0.219 


1.77 X 10^ 


2.36 


+1 


15.9 


0.253 


423 


0.218 


1.77 X 10^ 


2.33 


rUg = 200 GcV 


15.8 


0.248 


423 


0.217 


1.77 X 10^ 


2.34 


500 GoV 


15.9 


0.252 


423 


0.218 


1.77 X 10^ 


2.34 


1000 GcV 


15.9 


0.255 


423 


0.219 


1.77 X 10-^ 


2.34 


rUg = 200 GeV 


15.8 


0.251 


421 


0.214 


1.76 X 10^ 


2.29 


500 GeV 


15.9 


0.253 


423 


0.217 


1.77 X 10^ 


2.33 


1000 GeV 


15.9 


0.254 


424 


0.219 


1.77 X 10^ 


2.34 



Table 6: The NLO+NLL cross sections for stop-antistop production at the Tevatron and the LHC. 

Wc compare the SPSla' default input for the stop mixing angle and the light-flavour squark and 
gluino masses with various other choices for these SUSY parameters. Note that only one parameter 
is changed at a time, while the others are kept at their default values. The MSTW pdfs have been 
adopted and the scale has been set to the final-state stop mass. 

Cross sections for 6i6i production are collected in Table 7. We compare the LO bottom- 
quark-induccd contributions with the LO and NLO-I-NLL predictions based on the stop-like 
contributions that exclude bottom-quark initial states. Using the notation introduced in 

Sect. 2 and mi = m? — m? , the LO bottom-quark- induced contributions read 

y 01 

/ m?s [1 - cos(40^)] s + 2m?^ l + cos(40g) s + 2ml\ 
V 8{mjs + mi) 3i 4 ^ 2NcS )^ 

/mi + smj (s + 2m'i)[3 + cosi4e-^)]\ / 1 - P + 2ml/s \ 
^ V Ncs^ 8's ) ^^\l + P + 2ml/s) 

As these contributions depend on the gluino mass, we give results for rrig = 200, 500 GcV 
and 1 TeV. From the numbers presented in Table 7 it is clear that the bb — ^ bibi channel is 
always strongly suppressed, with cross sections well below 1% of the stop-like contributions. 



(0) _ asT^Cp 

bb-^bj>i ~ NrS 
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Tevatron (1.96 TeV) 


LHC @ 7 TeV 


LHC i 


I 14 TeV 




mi, [GeV] 


100 


200 


100 


400 


100 


400 


SPSla' default 

ATT 1 ATT T 


15.9 


0.253 




423 


0.218 




1.77 X 10 


2.34 




LO 


12.6 


0.227 




305 


0.156 




1.35 X 10 


1.67 




LO 56-clianncl only 


0.404 X 10-2 


0.330 X 10 


-4 


0.275 


0.346 X 10 


-3 


1.40 


0.564 X 10-- 


2 


LO 56-channcl only 
with mg = 200 GcV 


0.986 X 10^2 


0.870 X 10 


-4 


0.659 


0.667 X 10 


-3 


3.35 


0.111 X 10" 


1 


500 GcV 


0.454 X 10^2 


0.399 X 10 


-4 


0.309 


0.408 X 10 


-3 


1.58 


0.665 X 10" 


2 


1000 GeV 


0.335 X 10-2 


0.220 X 10 


-4 


0.227 


0.220 X 10" 


-3 


1.16 


0.360 X 10- 


2 



Table 7 : The LO and NLO+NLL cross sections for sbottom-antisbottom production at the Teva- 
tron and the LHC. We compare the default SPSla' prediction for the stop-like contributions with 
the LO contributions induced by bottom-quark initial states. The MSTW pdfs have been adopted 
and the scale has been set to the final-state sbottom mass. 



Bottom-quark-induced t-channel gluino exchange also leads to hih\ and h\b\ final states. 
The LO cross section for these processes is given by 



NrS 



1 



Nr_ 2m|s — mi -t- (2m|s -I- mi) cos(46'^ 



8N, 



[l-cos(40g)] + 



(mi -|- m|s) 



/3 



+ 



mi(l - cos(40(;)) + 4sm| (s + 2m'i){l - cos(4^^)) 



ANcs{s + 2mi) 



8s 



log 



1 - l3 + 2mi/s 
1 + l3 + 2m?_/s . 



with the identical expression for the charge conjugate process hh — > The corresponding 
numerical results for b-J)i production are listed in Table 8. Also the 6^6^ and h}J)k processes 
are suppressed by the small bottom-quark pdfs and never exceed the per-mille level with 
respect to 6^6^ production. 
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a{pp/pp -^hhi) [pb] 





Tevatron (1.96 TeV) 


LHC (§ 


7 TeV 




LHC @ 14 TeV 


mi, [GeV] 


100 


200 


100 


400 


100 


400 


SPSla' default 


0.111 X 10-2 


0.188 X 10- 


-4 


0.716 X 10-1 


0.205 X 10 


-3 


0.362 


0.306 X 10 


-2 


m§ = 200 GeV 


0.568 X 10-2 


0.518 X 10- 


-4 


0.335 


0.242 X 10 


-3 


1.64 


0.376 X 10 


-2 


500 GeV 


0.157 X 10-2 


0.247 X 10- 


-4 


0.994 X 10-1 


0.234 X 10 


-3 


0.500 


0.349 X 10 


-2 


1000 GeV 


0.447 X 10-3 


0.846 X 10- 


-5 


0.297 X 10-1 


0.124 X 10 


-3 


0.153 


0.187 X 10 


-2 



Table 8: The LO cross sections for bibi production at the Tevatron and the LHC. The MSTW 
pdfs have been adopted and the scale has been set to the final-state sbottom mass. 
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